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We study the effects of interband hybridization within the framework of an extended Falicov- 
Kimball model with itinerant c and / electrons. An explicit interband hybridization breaks the 
U(l) symmetry associated with the conservation of the difference between the total number of 
particles in each band. As a result, the degeneracy between multipolar electric and chiral orderings 
is lifted. We analyze the weak- and strong-coupling limits of the c-f electron Coulomb interaction at 
zero temperature, and derive the corresponding mean-field quantum phase diagrams at half-filling 
for a model defined on a square lattice. 
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I. INTRODUCTION 

The Falicov-Kimball modei^ (FKM) was primarily in- 
troduced to describe the metal-insulator transition of the 
mixed- valence compound SmBg. Later on, the model 
became widely accepted minimal Hamiltonian for 
studying several strongly correlated electron systems^ - — 
in particular, heavy fermion compounds^— In its origi- 
nal form, the FKM contains an itinerant c band of elec- 
trons that interact via a local Coulomb repulsion with 
localized / electrons. The spin degree of freedom of the 
electrons is not included. The local / electron number is 
strictly conserved and c-f electron coherence cannot be 
established^ An explicit hybridization between / and c 
orbitals provides an opportunity to overcome this short- 
coming i 13 ' 14 More recently, it was shown that a finite / 
electron bandwidth also induces c-f electron coherence, 
i.e., it can lead to an excitonic condensate even in absence 
of an explicit interband hybridization ! 15 ' 16 

These extended versions of the FKM were used to sub- 
stantiate the exciting idea of electronic ferroelectricity 
(EFE)f^— The ferroelectric phase only appears when 
the c and / orbitals have opposite parity under spatial 
inversion. The concomitant spontaneous breaking of in- 
version symmetry results from a nonvanishing average of 
(c'f). Since this expectation value corresponds to (exci- 
tonic) pairing of electrons and holes from different bands, 
the appearance of EFE is directly related with the for- 
mation of an excitonic insulator (EI) .— ~— 

The FKM with two dispersive bands, the so-called 
extended Falicov-Kimball model (EFKM), was studied 
previously for describing different properties of the EI 
phasej22r— However, as it was shown for the case of 
opposite-parity orbitals ; 15 ' 16 the inclusion of a finite in- 
terband hybridization can be very relevant because it re- 
moves the U(l) symmetry associated with the conserva- 
tion of the difference between the total number of parti- 
cles in each band: N c —Nf. In particular, this hybridiza- 



tion term is certainly relevant when the ground-state of 
the EFKM corresponds to an excitonic condensate. For 
the EFKM with interband hybridization (HEFKM), the 
excitonic condensate is in general replaced by Ising-type 
phases that only break discrete symmetries of the Hamil- 
tonian. 

In this paper we present a mean-field study of the influ- 
ence of an explicit hybridization on the symmetry-broken 
states that can take place for the HEFKM. To determine 
the ground-state quantum phase diagram of the HEFKM 
in the strong- and weak-coupling limits of c-f electron 
interaction, we assume that the interband hybridization 
amplitudes are small compared to the intraband hopping 
(transfer) integrals. Given the nature of the discrete sym- 
metries of the HEFKM, the natural ground-state candi- 
dates are chiral phases (CHPs) and states with multipo- 
lar electric orderings. 



II. MODEL 

By expressing the orbital flavor as a pscudospin vari- 
able, c\ = cL and fj = cj^ , the Hamiltonian takes the 
form 



(ij), 



u n A n a + v X! ( c !t c 4 + H 
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Here (ij) indicates that i and j are nearest-neighbor 

sites. The fermionic operators c^J annihilate (create) 
an electron on the spin a Wannicr orbital of the lattice 
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site Ri. The lattice has a total number of N sites, and 
nj a = c ]cr c ja i s the particle number operator for site j 
(a = {t;l})- £ a denotes the on-site energy for each or- 
bital, t a are the intraband hopping amplitudes, U is the 
local interorbital Coulomb interaction strength, and i> 7 
are the interband hybridization amplitudes, where 7 = 
for on-site hybridization and 7 = {tl, It} f° r intcrsite 
hybridization. The EFKM is recovered from Eq. (JTJ) by 
setting v-y = 0. In this limit, the model has a continuous 
U(l) symmetry, which is removed by the inclusion of an 
explicit hybridization. The discrete symmetries that re- 
main for the more general HEFKM are spatial inversion 
and time-reversal invariance. 

The pseudospin language of Eq. ([T]) unveils the sim- 
ilarity of H with other generic many-body Hamiltoni- 
ans. The EFKM (v 1 — 0) becomes an asymmetric Hub- 
bard model, i.e., a single band model for electrons with 
a spin-dependent dispersion. We will still use the name 
" Falicov-Kimball model" to indicate that the pseudospin 
degree of freedom represents a physical orbital degree of 
freedom. From now on, we will consider that T-L is de- 
fined on a square lattice and (n^ + tin) = 1 (half-filled 
band case) . We will also restrict to zero temperature and 
measure all energies in units of = 1. Finally, we will 
assume that the Wannier functions of c and / orbitals, 
</^(r — Rj) and ^.(r — Rj), are real. 

III. ORDER PARAMETERS 

In the rest of the paper we will refer to the pseudospin 
simply as "spin." The spin representation used in Eq. (TTJ) 
unveils the SU(2) structure of this internal degree of free- 
dom. This degree of freedom is the only one that survives 
at low energies in the large U/\t a \ limit. Consequently, 
the three different local or real-space order parameters 
correspond to the three components of the local spin vari- 
able, 

S i = ^Ei^'V' ( 2 ) 

a. a' 

where <r is the vector of the Pauli matrices. More com- 
plicated (or higher order) real-space order parameters in- 
volve products of spin operators in more than one unit 
cell. 

A real space modulation of (<SJ) leads to orbital order- 
ing. Here we will only consider the ordering wave vector 
Q = (71", 7r) that leads to staggered orbital ordering (SOO) 
because the effective interaction is antifcrromagnctic be- 
tween nearest-neighbors and the lattice under consider- 
ation is bipartite. The corresponding order parameter 
is 

^soo = E eiQ ' Rj ^) ■ ( 3 ) 

3 

If the two orbitals have opposite parity, a nonzero (ST) 
implies the presence of a spontaneous local electric po- 
larization that turns out to be uniform for the HEFKM. 



This is the EFE that was found in previous works for par- 
ticular limits of the HEFKM^ - — The uniform electric 
polarization is given by 

(V)= P J2(S?) (4) 

3 

with the interband dipole matrix element 

p = 2e J d 3 r t (v) r ^(r) , (5) 

where e is the electron charge. This phase breaks the 
spatial inversion symmetry of %. 

If the two orbitals have the same parity (for instance 
s and d orbitals), a nonzero modulation of (<SJ) corre- 
sponds to an electric quadrupole density wave (EQDW) 
as long as the tensor 

Ca> =e[ d 3 r</v(r)rv7v<Mr) (6) 

is nonzero for a' = a = — a {v,v' = {x,y, z}). In case 
the tensor q CT5 . [Eq. ©] vanishes, one has to look for 
the lowest order electric multipole that has a nonzero 
matrix element between the orbitals </>f(r) and </>j.(r). In 
the rest of this paper, we will assume that the tensor 
q<j ff does not vanish. In second quantization, the local 
electric quadrupole tensor on the unit cell j is given by 
the operator, 

Qj = qtt n Jt + qu^-i + 2 qnSj . (7) 

We note that q^ = q^. The corresponding quadrupolar 
order parameter in momentum space is given by 

3 

Again, for the Hamiltonian under consideration, the wave 
vector of the electric quadrupolar ordering is k = Q. 
Equation implies that a nonzero modulation of the x- 
spin component, (S®), corresponds to an EQDW. While 
translational symmetry is broken in this phase, time- 
reversal and spatial inversion symmetries are conserved. 
We note that the first two terms of Eq. (J7J) imply that 
orbital ordering will also lead to an EQDW. However, 
as it is also clear from Eq. (J7|), the quadrupolar ten- 
sors associated with the ordering along the z and x axes 
are different. The quadrupolar electric moment that is 
modulated under staggered orbital ordering (staggered z 
component) corresponds to a linear combination of the 
tensors q-|~|- and q^. On the other hand, the staggered 
ordering of the a;-spin component involves a modula- 
tion of a quadrupolar electric tensor proportional to q^ 
(hybridization- induced quadrupolar electric moment). In 
order to simplify the notation, we will use "EQDW" to 
denote the staggered ordering of the x-spin component 
and SOO for the the staggered ordering of the z-spin 
component. 
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Finally, a nonzero (<Sj) implies the spontaneous emer- 
gence of a current-density distribution between the two 
orbitals of the unit cell j . 21 ' 27 This can be easily verified 
if we neglect the overlap between orbitals that belong 
to different unit cells. In that case, the current-density 
operator at a point r near Rj is given by 

j,(r) = AsJ^^r-R^Vr^r-R,-), (9) 

e a 

where m e is the electron mass and the prefactor a takes 
the value +1 for j" and — 1 for 4-. We note that this current 
density flows between the two orbitals of the same unit 
cell (these are atomic currents when the two orbitals be- 
long to the same ion), in contrast to the orbital currents 
found in Ref. [28| that flow between different unit cells. 
This CHP breaks time-reversal symmetry and the physi- 
cal order parameter is the lowest order nonzero multipolc 
of the current-density distribution given by Eq. (j9]). The 
chiral ordering is staggered for orbitals with the same par- 
ity and uniform for orbitals with opposite parity. For in- 
stance, if we are considering two p orbitals, the staggered 
chiral ordering of the y-spin component corresponds to 
orbital antifcrromagnctism, because the current distri- 
bution given by Eq. §§§ generates a net magnetic dipolc 
moment. Since we will consider the general case of any 
arbitrary pair of orbitals, we will use "CHP" to denote 
the uniform ordering of the y-spin component (same par- 
ity orbitals) and "staggered chiral phase" (SCHP) to de- 
note the staggered ordering. 

IV. STRONG-COUPLING REGIME 

For the case | to- 1 , \ v y \ <C U, we can perform a large- U 
expansion, thereby reducing the HEFKM [Eq. ([T])] to an 
effective strong-coupling Hamiltonian, H sc , that repro- 
duces the low-energy spectrum of the original model. 

A large on-site Coulomb interaction splits the spec- 
trum of the HEFKM Hamiltonian into high- and low- 
energy parts. For t a = w 7 = 0, the lowest-energy sub- 
space is generated by the 2 N states that have one elec- 
tron per site. The high-energy subspaces are separated 
by energy gaps equal to Urid, where rid is the number of 
double occupied sites. For nonzero t a and u 7 , the elec- 
trons are no longer completely localized at their ions, i.e., 
an electron can gain kinetic energy by visiting virtually a 
neighboring site. Since we are considering the half-filled 
band case (one particle per site) , the low-energy effective 
model becomes a spin Hamiltonian T-L sc . The expression 
for T-L sc up to second order in the kinetic-energy terms is 

(ij) 

+ + J ZX S*S* + C) 

m 

+ £(BS i * + 2« b Sf), (10) 



where 
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Jxx =jj mh + v nm) 


(11) 


4 , 

J yy = 77 (*t*4. - V U V W 

U 


(12) 


2 

Jzz =77 (if + 1± — - 
u 


-«*t) > ( 13 ) 


4 , 

Jxz = jj (t t v lt - ti,v n ) 


(14) 


4 , 

Jzx = Jj (t t v n - tlVtf) 


(15) 


C=e t + £l -^j (t 2 + t\ + v 2 n + « 2 t ) , (16) 

B =e t - £ l • ( 17 ) 

It is well known that the half-filled isotropic Hubbard 
model can be mapped on an effective Hciscnbcrg model 



in the limit of a large Coulomb repulsion. For the more 
general EFKM, the intraband hopping amplitudes and 
the different on-site potentials lead to an effective XXZ 
model in a magnetic field B along the z axis. B is simply 
the energy difference between the two orbitalsJ^ As ex- 
pected, the interband hybridization of the HEFKM gen- 
erates anisotropic terms that explicitly break the U(l) 
invariance under uniform spin rotations about the z axis. 
While the intersite hybridization leads to anisotropic ex- 
change terms, the on-site hybridization leads to a Zeeman 
coupling to a uniform field along the x axis. 

For low enough values of B and no interband hybridiza- 
tion, the ground-state of H sc exhibits SOO. The simple 
reason is that J zz > \ J X x.yy\, i-e., the effective XXZ model 
is easy-axisJ^ If J zz is significantly larger than \ J X x.yy\, 
the SOO remains robust when the interband hybridiza- 
tion is included. Clearly, there exists a critical value of B 
that leads to a spin-flop transition to an ordered phase in 
the XY plane with a uniform component along the z axis 
(canted XY phase). In absence of interband hybridiza- 
tion, the U(l) invariance of the EFKM implies that spin 
component perpendicular to the applied field can point 
along any direction of the XY plane. In other words, 
there is a continuous ground-state degeneracy that in- 
cludes the EQDW (EFE) and SCHP (CHP) for orbitals 
with the same (opposite) parity. In this case, the inclu- 
sion of interband hybridization is very relevant because 
it lifts the continuous degeneracy and stabilizes only one 
of the two possible Ising-type ordcrings (along the x or y 
spin direction). 

For orbitals with opposite parity, we have t-\t\. < 0, 
= —Vfi and vq = 0. These relationships are de- 
rived from simple symmetry considerations. In this case 
Jxx = -z?(\tfh\ + tfi), Jyy = -£(\tth\ - u ?4.)> and 
Jzx = - Jxz- Since J X x,J y y < and \ J XX \ > \ J yy \, the 
energy is minimized by a ferromagnetic alignment of the 
spins along the x direction that corresponds to an EFE 
phase. Since this was previously shown in Ref. HH from 
now on we will concentrate on the case of equal parity 
orbitals. In this case, we have tft± > 0, v^f = v^, and v 
can be nonzero if the two orbitals belong to different ions. 
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Then J xx > Jyy and J zx — J xz - Since J xx > Jyy > 0, 
Vf± favors staggered Ising-type ordering along the x di- 
rection, while vq favors a uniform polarization along the 
x direction and, consequently, a staggered Ising-type or- 
dering along the y direction. In other words, the inter- 
band hybridization can stabilize an EQDW or a SCHP 
depending on the ratio between the on-site and intcrsitc 
hybridization amplitudes. 

We introduce now the mean-field variational states and 
the corresponding energies for the three order parameters 
that we introduced in the previous section. 

(i) SOO. This phase has a staggered spin component 
along the z-direction, and a uniform component along 
the x direction that can be induced by the on-site hy- 
bridization term vo-, 



= S(sin0i,O,cos0ie* QR >") , (18) 
DNS 2 J XX sin 2 0i 

(19) 



E$°° = - DNS 2 J ZZ cos 2 6»i 4 
+2v SNsm6 1 + DNC 



(ii) SCHP. This phase has a staggered spin component 
along the y direction and uniform polarizations along the 
x and z directions, 

{Sj)e x e 2 = S (sin 6>i cos 2 , e lQR3 sin Q x sin 2 , cos X ) , 

(20) 

£ SCHP = DNS 2 ((J xx + Jyy) cos 2 2 (1 - cos 2 6i) - Jyy 
+ (Jyy + J zz ) cos 2 Q\ + 2J XZ sin 9\ cos 0i cos 6 2 ) 
+DNC + NSB cos 9-i + 2NSv sin 9 1 cos 9 2 . 

(21) 

(iii) EQDW. In this case the staggered spin component 
is aligned along the x direction and there is a uniform 
component along the z direction induced by B (the y 
component vanishes), 

(Sj) ei =S(e lQR > sin 0i,O, cos 0i) , (22) 
^eqdw = DNC + NBScosdl _ DNJ XX S 2 sin 2 0i 

+DNJ ZZ S 2 cos 2 9± . (23) 

In all cases we have D = 2, Q = (it, tt), and S = 1/2. 
By minimizing the respective energies with respect to 0i 
and 02, we determine the quantum phase diagram as a 
function of the band-structure parameters. 

Figure [Tla) shows the effect of a finite intersite hy- 
bridization amplitude, v^, on the quantum phase dia- 
gram of the EFKM. The intersite hybridization stabilizes 
the EQDW relative to the SCHP. On the other hand, the 
EQDW is also favored relative to the SOO because v^i 
decreases the value of J zz and simultaneously increases 
the value of J xx [see Eqs. (JTTJ) and ([T3"])]. Figure [T^b) 
illustrates the effect of a finite on-site hybridization. In 
this case, the SCHP is favored relative to the EQDW. In 
contrast to u-j-j., vq does not change the transition point 
between the SCHP and the SOO. The simple reason is 
that wo does not affect the exchange constants. It just 
generates an effective pseudomagnetic field along the x 
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FIG. 1: (Color online) Ground-state phase diagram of the 
2D EFKM in the strong-coupling regime. Band-structure pa- 
rameters are = 0.0, tf = 1.0, and U — 10. Left-hand side 
diagrams [in panels (a) and (b)] give results for the nonhy- 
bridized EFKM (vo = and u^j. = 0), while right-hand side 
diagrams show the dependence on the (a) intersite hybridiza- 
tion v-fi and (b) on-site hybridization vo for e-j- = 0.15. Panel 
(c) gives the stability region of the staggered chiral phase and 
the electric quadrupole density wave in dependence on vo and 
Vfi for e-f = 0.5, t± = 0.5. 



axis that leads to a finite canting angle in both phases. 
Figure [He) shows the phase diagram as a function of vo 
and v-fi for a large enough value of B = 0.5 and same 
parity orbitals. Again, we can see that a finite on-site 
hybridization vq strengthens the SCHP, while the inter- 
site hybridization v^i stabilizes the EQDW. Within our 
simple mean-field approximation the boundary between 
these two phases is a straight line. 
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WEAK-COUPLING REGIME 



with 



It has been shown in Ref. [l8| that the mean-field 
ground-state phase diagram of the 2D EFKM agrees al- 
most perfectly with the one obtained by a constrained 
path Monte Carlo technique, even in the intermediate 
coupling regime. This agreement motivates us to perform 
a Hartrcc-Fock decoupling of the HEFKM to explore the 
quantum phase diagram for small U/\t a \. 

The weak-coupling analysis requires to express the rel- 
evant order parameters in momentum space. In partic- 
ular, the Fourier components of (<SJ) and (<SJ) can be 
represented as a complex number, 



i^qi^Se^ 



N N--k+Qt L ki/ 
k 



u 



u 

N 



where 



3 

Vn 4" 3 



(24) 
(25) 
(26) 



(27) 
(28) 



£kcr =£cr + Un-o + 2i cr (cosfc :E + cos k y ) , (30) 
«k =vo + 2u t 4.(cos k x + cos k y ) , (31) 



TV 

k 



A* - 
Q _ 



kJ, C k+Qt/ 



(32) 
(33) 
(34) 



The mean-field Hamiltonian (|29| can be easily diagonal- 
izcd by the canonical transformation^ 

Ck.m = u k,rnC k -|- + Uk,m c k4. ^ k > mC k+Qt- ~l~ ^k,mC k _^Q^ , 

(35) 

where m = 1, 2, 3, 4. The coefficients arc solutions of the 
associated Bogoliubov de Gennes equations, H^^^m = 
^k,m*k,m , with 



1/WC _ 



/ e kt « k -AJ, \ 

"k £ki -A Q 

-A^ e k+Qt v k+Q 

\-A Q w k+Q e k +cu/ 



(36) 



and * kim = (ii k ,m ! Wk,m ! '"k,TO ) «k,m)' r - The energy per 
site results as 



pwc 
^0 



N N 



=ij2' E ^f(EKm) - Un tni + i|A Q | 2 , (37) 



The ordering wave vector Q determines the modulation 
of the real-space order parameter. According to our 
strong-coupling analysis, we have Q = (tt, it) for orbitals 
with the same parity and Q = (0, 0) for orbitals with op- 
posite parity. Again, the U(l) invariancc of the EFKM 
(u 7 = 0) implies that the energy does not depend on 
<p. Consequently, there is an infinite number of ground- 
states with Aq 7^ that results from the spontaneous 
U(l) symmetry breaking of the EFKM (excitonic conden- 
sate). A finite intcrband hybridization (i> 7 ^ 0) removes 
the continuous U(l) symmetry and lifts the tp degeneracy 
of the HEFKM ground-state. 

For orbitals with opposite parity, the hybridization in 
momentum space takes the form i; k = 2iv^^(suik x + 
sin /c y ) and the EFE state has a lower energy than the 
CHP. This is in agreement with the result for the FKM 
extended by a (small) intersite hybridization in Ref. Il4l . 
Therefore, from now on we will focus only on the equal 
parity case. The Hartree-Fock decoupling suggested by 
Eq. ([M]) gives 



H wc =^e kCT cLc k(T 



k<r c k- 



Q C k;i c k+Qt 



E a q4 



k+Qt L ki 



(29) 



where /(-E kjm ) is the Fermi function containing the new 
quasiparticle energies -E k , m and the prime denotes that 
the k summation extends over the magnetic Brillouin 
zone only. The chemical potential fi is determined by 
the condition 



^EW 



k .lit J 



(38) 



k,m 



Next we consider the mean-field decoupling that leads 
to SOO. In this case, we introduce the possibility of a 
periodic modulation in the electronic density with inde- 
pendent amplitudes for each spin polarization, 



8 a cos(QR 4 ), 



with 



(39) 



(40) 



The associated Bogoliubov de Gennes equations are 

nl°°^ Km = E^o^ k m: with 



t/SOO 

rt k 



/£kt Vk USf 

Vk Ski u h 

US-f £k+Qt v k+Q 

\ USl Wk+Q £k+Ql/ 



(41) 
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The SOO order parameter becomes 



^soo 



(42) 



For asymmetric bands, t-j- 7^ tj,, the presence of a nonzero 
SOO leads to a secondary charge-density-wave (CDW) 
order, whose order parameter is given by 



Jcdw 



5f + 5i 



(43) 



This secondary CDW provides a simple way of detecting 
the SOO in real materials. The mean-field energy per 
site that results from such a kind of SOO is 



N 



1 Y^'p 



SOO f( pSOON 



k,m 



(44) 



By solving the self-consistency Eqs. (|33| and (j40|) . and 
comparing the corresponding mean-field energies given 
by Eqs. (|37| and (|44|) . we compute the ground-state 
phase diagram. 

Figure [2] is the weak-coupling counterpart of Figure [T] 
Like for the strong-coupling regime, Fig. [5{a) shows that 
an increasing value of narrows the SOO phase while 
the region of the EQDW phase is enlarged. On the other 
hand, the on-site hybridization vq favors the SCHP rel- 
ative to the EQDW and it does not have a noticeable 
effect on the transition line between the SOO and the 
SCHP [see Fig.^b)]. This also coincides with the strong- 
coupling results. Figure [2jc) shows the stability regions 
of the EQDW and SCHP as a function of the hybridiza- 
tion amplitudes for a large enough — ejj = 0.5. Quali- 
tatively, the result of our Hartrec-Fock approach is simi- 
lar to the one obtained from the strong-coupling analysis 
[see Fig. [He)]. However, a more quantitative analysis 
shows that the area of stability for the SCHP is reduced 
relative to the strong-coupling result. A large Coulomb 
repulsion inhibits hopping processes and consequently re- 
duces the influence of the intersite hybridization rel- 
ative to the effect of the on-site hybridization Vq. 
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FIG. 2: (Color online) Ground-state phase diagram of the 
2D EFKM in the weak-coupling regime. Band-structure pa- 
rameters are e± — 0.0, t-f = 1.0, and U = 2. Left-hand side 
diagrams [in panels (a) and (b)] give results for the nonhy- 
bridized EFKM (vo = and Vf± = 0), while right-hand side 
diagrams show the dependence on the (a) intersite hybridiza- 
tion Vfi and (b) on-site hybridization vo for = 0.15. Panel 
(c) gives the stability region of the staggered chiral phase and 
the electric quadrupole density wave in dependence on vo and 
Vfi for e-f = 0.5, t± = 0.5. 



VI. CONCLUSIONS 



The Hamiltonian that we considered in this work is a 
very simple extension of the Falicov-Kimball model. In 
spite of its simplicity, we have shown that this model 
leads to a very rich quantum phase diagram that con- 
tains all the possible local order parameters (three dif- 
ferent components of the local spin Sj) considered in 
Sec. III. The ordering wave vector Q is selected by the 
nesting property of the noninteracting Fermi surface in 
the weak-coupling limit and by the antiferromagnetic na- 
ture of the exchange interactions on a bipartite lattice in 
the strong-coupling limit. Most notably, the stability of 
the different broken symmetry states is very sensitive to 



a few band-structure parameters. According to these re- 
sults, it is necessary to have very accurate information 
about the band-structure properties near the Fermi en- 
ergy to predict the correct ordered state. In particular, if 
the two orbitals have different angular momentum (like s 
and d orbitals), the SCHP may remain hidden to most of 
the experimental probes. The simple reason is that the 
spontaneous current-density distribution given by Eq. © 
has no net magnetic moment. Consequently, this phase 
can only be detected by using an experimental probe that 
couples to the lowest nonzero multipole of the current- 
density distribution. The SCHP becomes stable above 
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a critical value of the on-site hybridization as long as 
the diagonal energy difference between the two orbitals 
\s c — £f \ is also larger than a minimal value. 

Although all the calculations of this work were done 
for D = 2, we do not expect any qualitative change for 
D > 2. The obtained consistency between the weak- 
and the strong-coupling approaches suggests that our 
results are robust. In particular, the absence of geo- 
metric frustration in the strong-coupling regime, whose 
weak-coupling counterpart is the nesting property of the 
Fermi surface, facilitates the search for the broken sym- 
metry state that minimizes the energy for each set of 
Hamiltonian parameters. For orbitals with opposite par- 
ity under spatial inversion, we confirmed that the ferro- 
electric phase has always a lower energy than the chi- 
ral phase. For orbitals with the same parity, we found 
that the stabilization of the electric quadrupole density 



wave or the staggered chiral phase depends strongly on 
the dominant interband hybridization. The on-site hy- 
bridization, that is only allowed when the two orbitals 
belong to different ions, favors the staggered chiral phase, 
while a nearest-neighbor interband hybridization favors 
the electric quadrupole density wave. 
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